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Abstract 

We study the spectrum of confining strings in SU(3) pure gauge 
theory, by means of lattice Monte Carlo simulations, using torelon 
operators in different representations of the gauge group. Our results 
provide direct evidence that the string spectrum is according to pre- 
dictions based on n-ality. Torelon correlations in the rank-2 symmetric 
channel appear to be well reproduced by a two-exponential picture, 
in which the lowest state is given by the fundamental string ai = a, 
the heavier string state is such that the ratio 02! (y\ is approximately 
given by the Casimir ratio Cgym/Q = 5/2, and the torelon has a much 
smaller overlap with the lighter fundamental string state. 



The spectrum of confining strings in 4-d SU(A^) gauge theories has been 
much investigated recently. Several numerical studies in the context of a 
lattice formulation of the theory have appeared in the literature, providing 
results for color sources associated with representations higher than the fun- 
damental, see e.g. Refs. [H I2l El 13 El E] • General arguments show that the 
string tension must depend only on the n-ality, k = mod{l,N), of a repre- 
sentation built out of the (anti-) symmetrized tensor product of / copies of 
the fundamental representation. The confining string with n-ality k is usu- 
ally called fc-string, and ak is the corresponding string tension. Using charge 
conjugation, ak = CN-k- As a consequence, SU{3) has only one independent 
string tension determining the large distance behavior of the potential for 
k ^ 0. One must consider larger values of to look for distinct /c-strings. 
Lattice results for = 4, 5, 6 [HJ E] show a nontrivial spectrum for the k- 
strings. In particular the data of Ref . [H] , obtained using color sources in the 
antisymmetric representations of rank k, turn out to be well reproduced by 
the sine formula 

Ofc ^ sm{kiT/N) 

a ^ sin(7r/Ar) ' ^ 

[a = ai) within their errors. The sine formula has been suggested by several 
theoretical works, especially in the context of supersymmetric theories, see 
e.g. Refs. 13 E] and references therein. 

On the other hand, numerical results for different representations with 
the same n-ality apparently contradict the picture that n-ality is what really 
matters. For example, in the SU(3) case, Monte Carlo data for the Wilson 
loops for several representations |21 E] show apparently area laws up to rather 
large distances, approximately 1 fm, also for representations with zero n-ality, 
and the extracted string tensions turn out to be consistent with the so-called 
Casimir scaling |I] . In the lattice study of Ref. E] ? considering larger val- 
ues of N, the fc-string tensions were extracted from the torelon masses, i.e. 
from the exponential decay of correlations of characters of Polyakov lines. 
In Ref. jH], while the antisymmetric representations provided rather clean 
measurements of ak reproducing the sine formula, the numerical results for 
the symmetric representations suggested different values of the correspond- 
ing string tensions. For example, in the case of rank 2, o"sym/c" ^ 2, which is 
approximately the value suggested by Casimir scaling or by the propagation 
of two noninteracting fundamental strings. These results that apparently 
contradict n-ality have been recently discussed in Ref. 0. They have been 
explained by arguing that standard color sources, such as Wilson loops and 
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Polyakov lines, associated with representations different from the antisym- 
metric ones have very small overlap with the stable fc-string states, being 
suppressed by powers of in the large- limit, and in some cases also 

exponentially. Since A'" = 3 is supposed to be already large, these arguments 
may explain why the predictions of n-ality have not been directly observed 
in the numerical simulations, which are limited in accuracy. Moreover, the 
situation worsens in the case of larger N . 

Motivated by this recent work, we decided to return on this issue in the 
context of the 4^d SU(3) gauge theory. Performing Monte Carlo simulations 
of the SU(3) gauge theory in its Wilson lattice formulation, we measure 
correlators of Polyakov lines in the representations of rank 1 (fundamental) 
and 2 of the SU(3) group, in order to check whether their large- distance 
behaviors, and therefore the values of the corresponding string tensions, are 
consistent with ra-ality. 

In our numerical study we use a method based on torelon correlators PU] ■ 
The string tensions are extracted from the large-time behavior of "wall-wall" 
correlators of Polyakov loops in spatial directions, closed through periodic 
boundary conditions (see e.g. Refs.[TUll^): 

Gr{t) = E (Xr[P(0,0;0)] Xr[P{xi,X2-M. (2) 

where P(xi,a;2;t) = IixzU3{xi,X2,x^;t). f/(x; t) are the usual link variables, 
and Xr is the character of the representation r: Xi[P] = TrP for the funda- 
mental representation, while the two representations of rank 2 (antisymmetric 
and symmetric, both with n-ahty k = 2) have Xasym[-P] = | ((Tr P)^ — Tr P^), 
Xsym[P] = |((TrP)2 + TrP2), respectively. Note that, since for SU(3) 
Xasym[^] = Xfff^]) the Correlators in the k = 1 fundamental and k = 2 
antisymmetric representations are identical. We have also studied the ad- 
joint representation (n-ality k = 0), for which: Xa,di 

[P] = |TrPp - 1. In this 
case, one should consider the connected correlator, since (xadj[P]) 7^ 0. 

The correlators ^ decay exponentially as exp(— mfct), where ruk is the 
mass of the lightest state in the corresponding representation. Actually, on 
a finite lattice with periodic boundary conditions we have Gr{t) oc cosh(t — 
T/2), whete T is the temporal size. For a k-\oop of size L, the A;-string 
tension is obtained using the relation [TU] 

ruk = CTkL - (3) 
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The last term in Eq. is conjectured to be a universal correction, and it is 
related to the universal critical behavior of the flux excitations described by 
a free bosonic string [llj. 

We present results obtained at j3 = 5.9 and for two asymmetric lattices 
12^ X 24 and 16^ x 24, allowing us to compare the results using Polyakov lines 
with different length, i.e. L = 12, 16. Using our data for the fundamental 
string tension, see below, and the standard value y/a = 440 MeV, L = 12, 16 
correspond to approximately 1.5 and 2 fm, respectively. In our simulations we 
upgraded the SU(3) variables by alternating microcanonical over-relaxation 
and heat bath steps, typically in a 4:1 ratio. More details on the algorithm 
can be found in Ref. 0. We collected rather high statistics, ~ 16M sweeps 
(considering a sweep as the upgrading of all links of the lattice independently 
of the algorithm) for L = 12 and ~ 7M sweeps for L = 16. Measurements 
were taken every 20 sweeps. In order to improve the efficiency of the mea- 
surements we used smearing and blocking procedures (see e.g. Refs. P^) to 
construct new operators with a better overlap with the lightest string state. 
We constructed new super-links using three smearing, and a few blocking 
steps, according to the value of L, i.e. two for L = 12 and four for L = 16. 
These super-links were used to compute improved Polyakov lines. Our im- 
plementation of smearing and blocking is as follows |6j: Smearing replaces 
every spatial link on the lattice according to: 



where V indicates the projection onto SU{N) and the sum only runs on 
spatial directions. Similarly, blocking replaces each spatial link with a super- 
link of length 2a: 



Ukix)^v\Ukix)Ukix + k) + af Ujix)Ukix + j)Ukix + j + k))Ujix + 2k) 



The blocking procedure can then be iterated n times to produce super-links 
of length 2"a. The coefficients as and aj can be adjusted to optimise the 
efficiency of the procedure. We constructed new super-links using as = aj = 



In Figs, m and |21 we show the wall-wall correlators as a function of the 
distance t in the cases of fundamental and symmetric representations, from 
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Figure 1: Correlator data for the fundamental (circles) and rank-2 symmetric 
(squares) representations, with L = 16, as a function of the temporal distance 
t. The top solid line is a hyperbolic cosine fit to data with t > 3, leading to 
an estimate of mi ; the bottom solid line is a hyperbolic cosine with the same 
exponent, which is well supported by the symmetric representation data for 
t > 3. The dashed line represents a fit using the two-exponent Ansatz 
Data from the rank-2 symmetric representation at one less blocking step is 
shown in diamonds. 

the runs with L = 16 and L = 12 respectively. The data for the correla- 
tor in the fundamental representation allow us to accurately determine the 
fundamental string tension, and the two lattices provide consistent results 
using Eq. (jH)), i.e. a = 0.0664(5) and a = 0.0668(3) respectively from the 
L = 16 and L = 12 runs (obtained by fitting results starting from distances 
t = 3, 4, respectively). On the other hand, such an agreement is not observed 
in the case of the symmetric representation. However, the L = 16 data for 
the symmetric correlator shows a clear evidence that its asymptotic behavior 
is controlled by the fundamental string; indeed, fitting the data for t > 3 we 
obtain cXsym = 0.070(4), in agreement with rz-ality. Although data at small t, 
t < 3, show a clear contamination by heavier states, in the symmetric rep- 
resentation case the overlap with the fundamental string state of the source 
operator, obtained by performing four blocking steps after smearing, appears 
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Figure 2: Correlator data for the fundamental (circles) and symmetric 
(squares) representations, with L = 12, as a function of the temporal distance 
t. Solid lines are exponential fits to data with t > 3, leading to markedly 
different values of . 

to be sufficient to show the actual asymptotic behavior before the signal dis- 
appears within the error. This is not observed in the L = 16 data using the 
source operator with three blocking steps (one less) and in the L = 12 data 
(where two blocking steps are employed). Up to the distances that we can 
observe before the signals die off into the noise, the correlators appear to be 
dominated by the propagation of a much heavier state, which would suggest 
Csym ~ 0.16, whose corresponding ratio cTsym/o' ~ 2.4 is rather close to the 
Casimir ratio of the two representations, i.e. 5/2. 

A simple interpretation of the behaviour of torelon correlations is provided 
by a picture based on two propagating states; indeed, the numerical results 
suggest: 

Gsym(t) ^ Cie-'"^* + 026-^^' (6) 

where rrii = cTiL — 7r/(3L). Due to n-ality, ai = a. The string tension of the 
first excited string state is (72, and the overlap with the states satisfy ci -C C2. 
The data on the smaller lattice L = 12 already suggests that the ratio a2/(J 
should be approximately equal to the Casimir ratio. As shown in Fig. ^ the 
Ansatz (jH)) fits well the L = 16 data at four blocking steps for t > 1, with 
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oijo — 2.2 and Ci/c2 — 0.12. Assuming a mild dependence of q on L, and 
given the smallness of Ci , the signal from the fundamental string should only 
be visible for large enough t: 



Clearly, for an acceptable signal-to-noise ratio, t must not grow excessively, 
and the above equation requires a sufficiently large spatial size L ; this is 
consistent with the results from the lattice sizes we use. 

According to the large- arguments of Ref. P, the ratio Ci/c2 should be 
suppressed at least by a power of l/N"^. This would make the observation of 
the asymptotic k = 2 string state, using sources in the symmetric represen- 
tation, much harder for larger A^, thereby explaining the results for A^ = 4, 6 



In conclusion, the results that we have presented provide direct evidence 
that the spectrum of confining strings is according to predictions based on 
n-ality. Torelon correlations in the rank-2 symmetric channel appear to be 
well reproduced by a two-exponential picture, in which the lowest state is 
given by the fundamental string cxi = a, the heavier string state is such that 
the ratio o"2/o"i is approximately given by the Casimir ratio Csym/Cf = 5/2, 
and the torelon has a much smaller overlap with the lighter fundamental 
string state. 

In the case of the adjoint representation, since its n-ality is zero, the cor- 
responding mass of the exponential decay in the connected correlator should 
not depend on L, but it should be related to the propagation of gluelumps. 
Actually, since the gluelumps have a limited physical size, by increasing L 
we expect a smaller and smaller overlap of the adjoint Polyakov lines with 
the lowest states; this makes the evidence for the so-called adjoint string 
breaking^ very difficult when using the method applied here, requiring a pro- 
hibitive amount of statistics. Indeed, both the L = 16 and L = 12 data 
seem to correspond to the propagation of a string state with cXadj ~ 0.145, 
and therefore cXadj/o" ~ 2.1, which is again rather close to the corresponding 
Casimir ratio 9/4. 
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